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-We  introduce  general  framework  for  developing  and  evaluating  deter¬ 
ministic  models  of  production  processes.  We  evaluate  familiar  linear  pro¬ 
gramming  planning  formulations,  Manufacturing  Resources  Planning  (MRP) 
and  Critical  Path  Methods  (CPM)  in  terms  of  the  framework.  Flaws  are 
exposed  in  linear  programming  formulations  incorporating  lead  times.  We 
— provide  a  correct  reformulation,  as  well  as  a  more  general  formulation  accom¬ 
modating  non-integer  lead  times  and  unequal-length  planning  periods. 
Redefining  MRP  and  CPM  illustrates  the  framework's  use  for  comparing 
seemingly  unrelated  models  and  for  rigorously  explaining  limitations  in  exist¬ 
ing  models.  We"propose' a  structured  approach  for  the  development  of  accu¬ 


rate  models  of  production. 
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1.  INTRODUCTION 


Researchers  in  production  planning  and  scheduling  are  beginning  to  realize  that  their 
work  to  date  has  seen  a  disappointing  level  of  use  in  industrial  practice.  For  example,  Abra¬ 
ham  el  al.  [1985]  point  out  "a  disparity  between  the  insights,  orientation  and  methodologies 
of  past  [production]  research  and  the  possible  scientific  bases  of  current  practice."  They  stress 
the  need  for  a  "unifying  framework"  to  help  focus,  position  and  evaluate  future  research. 

Although  most  researchers  and  practitioners  would  agree  that  a  unifying  framework  is 
lacking,  there  is  little  consensus  as  to  what  it  should  consist  of  or  what  needs  it  should 
address.  Consider,  for  example,  the  following  proposed  frameworks.  Abraham  et  al.  [1985] 
suggest  a  hierarchical  framework  of  decision  levels  corresponding  to  a  presumed  hierarchy  of 
decision-making.  Dempster  el  al.  [1981]  propose  a  general  framework  for  the  analytical 
evaluation  of  hierarchical  planning  systems.  Bensoussan  et  al.  [1985]  offer  a  "unified 
mathematical  treatment  of  production  planning  and  production  smoothing  problems,  in  the 
framework  of  optimal  control  theory."  Akinc  and  Roodman  [1986]  propose  a  new  framework 
for  modeling  aggregate  production  planning  problems  in  which  emphasis  is  placed  on  flexibil¬ 
ity  in  specification  of  the  decision  options  and  the  relevant  cost  structure.  Geoffrion  [1985] 
proposes  a  "structured  modeling"  framework  to  facilitate  the  data-base  development  for  appli¬ 
cation  of  pre-specified  operations  research  models. 

Surely,  representation  of  the  appropriate  decision  alternatives,  decision  hierarchy,  and 
cost  structure,  as  well  as  the  effort  required  in  data-base  development,  are  important  concerns 
for  planning  and  scheduling.  However,  there  is  another,  equally  fundamental  concern  which 
attracts  our  research  efforts.  Embedded  in  each  planning  and  scheduling  model  is  a  model  of 
production  -  a  mathematical  representation  of  the  set  of  technologically  feasible  operations  or 
actions  within  the  production  process.  While  the  body  of  literature  proposing  optimal  or 
near-optimal  algorithms  for  production  planning  and  scheduling  models  is  voluminous,  a  very 
narrow  range  of  models  of  production  is  considered  in  this  literature,  as  well  as  in  the  above- 
cited  frameworks.  Clearly,  an  accurate  representation  of  the  production  process  is  critical  to 
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the  optimality  and  even  feasibility  of  generated  plans  and  schedules.  Therefore,  we  suggest 
that  a  framework  is  needed  to  (1)  assess  the  accuracy  or  validity  of  the  representation  of  a 
particular  production  process  by  a  given  production  model,  and  (2)  guide  the  development  of 
more  accurate  production  models  when  existing  models  provide  inadequate  representations  of 
the  production  process. 

In  this  paper,  we  introduce  a  general  framework  for  developing  and  evaluating  deter¬ 
ministic  models  of  production.  The  framework  incorporates  familiar  modeling  principles 
such  as  flow  conservation  and  activity  analysis,  but  it  expresses  tb:se  principles  in  terms  of 
new  abstract  model  elements.  When  existing  models  of  production  are  expressed  in  terms  of 
the  framework,  extensions  and  generalizations  become  apparent.  Moreover,  the  mathematical 
language  and  structure  of  the  framework  provides  a  foundation  to  compare  and  contrast 
seemingly  unrelated  models  arising  in  such  diverse  environments  as  discrete-parts  manufac¬ 
turing,  continuous-flow  production,  and  project  management. 

The  framework  is  introduced  in  Section  2.  In  the  sequel,  we  redefine  in  terms  of  the 
framework  three  of  the  most  familiar  production  planning  models:  standard  linear  program¬ 
ming  formulations.  Manufacturing  Resources  Planning  (MRP),  and  Critical  Path  Methods 
(CPM).  Each  model  highlights  different  abstract  model  elements  of  the  framework,  and 
serves  to  illustrate  different  uses  of  the  framework. 

In  section  3,  we  use  the  framework  to  develop  more  accurate  models  of  production. 
Familiar  linear  programming  (l.p.)  planning  models  are  analyzed,  including  those  which  incor¬ 
porate  production  lead  times".  We  consider  physical  phenomena  which  give  rise  to  lead 
time  and  show  that  the  physical  phenomena  should  be  categorized  into  three  distinct  types. 
Flaws  in  familiar  formulations  are  exposed,  and  we  derive  a  correct  reformulation  from  the 
framework.  We  then  generalize  the  formulation  to  accommodate  non-integer  lead  times  and 
unequal-length  planning  periods. 

Section  4  illustrates  use  of  the  framework  for  explaining  limitations  of  existing  models. 
We  show  that  MRP  techniques  are  based  on  an  incomplete  model  of  production.  As  a 


consequence,  a  simple  explanation  is  provided  of  why  MRP  techniques  will  in  general  lead  to 
excessive  work-in-process  inventories. 

Section  5  illustrates  the  value  of  a  coherent  language  for  comparing  production  model 
elements.  We  interpret  CPM  as  a  model  of  production.  Recognizing  that  a  strict  precedence 
relationship  between  two  activities  is  a  form  of  inventory  balance,  it  becomes  clear  how  CPM 
can  be  usefully  extended  by  incorporating  activity  analysis  of  l.p.  models. 

Section  6  concludes  the  paper  with  a  brief  outline  of  how  the  framework  can  be  used  in 
a  structured  modeling  approach  for  the  development  of  accurate  production  models. 


2.  A  GENERAL  FRAMEWORK 

Our  development  of  a  general  framework  for  models  of  production  follows  from 
research  it.  production  theory  initiated  by  R.W.  Shephard.  Shephard  [1970]  developed  the 
first  axiomatic  steady-state  model  of  production.  The  model  elements  in  his  framework  are 
production  correspondences  which  abstractly  model  the  input-output  possibilities  of  the  tech¬ 
nology.  The  first  general  activity  network  model  was  developed  in  Shephard  et  al.  [1977], 
which  introduces  intermediate  product  transfers  between  activities  and  activity  dynamic  pro¬ 
duction  correspondences.  Hackman  [1984]  extends  this  framework,  replacing  activity 
correspondences  with  activity  dynamic  production  functions.  Our  presentation  is  based  in 
part  on  Chapter  2  of  Hackman  [1984], 

The  framework  is  a  meta-mndr'  ction.  It  delineates  precisely  what  is  and  what 

is  not  a  valid  production  model  in  ,.iuch  the  same  way  that  the  "meta-model': 

min  t  v 

1  v  =  h 
x  >  0 

delineates  what  is  and  what  is  not  a  linear  program.  The  framework  identifies  the  abstract 
model  elements,  or  building  blocks,  and  the  relationships  between  *'  e  elements,  which  are 
required  in  a  mathematical  representation  of  a  production  system.  The  development  of  a 
deterministic  model  of  a  production  system  is  viewed  as  a  specification  of  the  elements  of  the 
framework  so  as  to  reflect  the  particular  technologies  and  processes  of  the  system.  We  now 
define  these  elements. 

The  Activity  Network 

The  framework  models  a  production  system  as  a  directed  network  whose  nodes 
represent  primitive  production  activities,  i.e. ,  activities  whose  internal  organization  is  not 
further  modeled.  Directed  arcs  indicate  possible  transfers  of  intermediate  products.  Inter¬ 
mediate  products  are  outputs  of  one  or  more  activities  which  serve  as  inputs  to  other  activi¬ 
ties.  Production  at  each  activity  requires  intermediate  products  transferred  from  other 


activities  and/or  system  exogenous  inputs.  System  exogenous  inputs  are  of  two  types:  non- 
storable  services  and  storable  materials.  Services  include  labor  trades,  machines  and  facilities, 
while  materials  include  purchased  parts,  raw  materials,  fuels,  etc.  An  exogenous  material  may 
be  the  same  as  an  intermediate  product  of  the  system. 

Allocations  of  system  exogenous  inputs  are  represented  as  transfers  from  an  initial  node 
,t0-  In  a  system  with  Ar  activities,  transfers  of  final  products  are  represented  as  transfers  to  a 
sink  node  .tv»i-  A  final  product  may  be  the  same  as  an  intermediate  product  of  the  system, 
as  in  the  case  of  spare  parts.  Thus,  a  production  system  is  regarded  as  a  jointly  operating, 

finite  number  of  interrelated  primitive  production  activities  2 . which  use  system 

exogenous  inputs  of  goods  and  services  to  produce  final  outputs. 

We  emphasize  that  the  activity  network  is  merely  an  abstract  representation  of  work 
flow  and  need  not  represent  any  physical  arrangement  of  facilities. 

Modeling  the  Flows  of  Goods  and  Services 

Since  flows  of  goods  and  services  (inputs  and  outputs)  have  a  dynamic  character,  we 
model  each  flow  as  a  bounded  function  of  time,  where  time  is  modeled  by  the  nonnegative 
part  of  the  real  line.  Time  0  is  a  reference  point  defining  the  point  in  time  after  which  flows 
would  be  determined  by  planning  calculations  using  the  model. 

There  are  two  fundamental  types  of  flows.  In  the  first  and  more  common  type,  called  a 
rate-based  flow,  x(t)  represents  the  rate--  quantity  per  unit  time--  of  flow  at  time  r.  In  the 
second  type,  called  an  event-based  flow,  x(t)  is  a  numerical  representation  of  an  event  at  time 
t.  For  example,  an  event-based  flow  in  a  project-oriented  production  system  might  be 

1  if  t  is  the  project  completion  time 
0  otherwise. 

In  another  example,  suitable  for  batch  transfers  of  intermediate  products.  ,v(r)  might  indicate 
the  quantity  transferred  at  time  r.  Unless  specifically  identified  as  event-based,  all  flows  in 
this  paper  are  rate-based. 


For  each  model,  an  index  set  A  =  {/*}f  of  epochs  is  defined.  The  /*’s  are  ordered  so 
that  0  =  to<t i<  .  Each  Ik=(tk-\-lk]  is  a  period,  the  points  in  A  are  called  time-grid 
points,  and  A  is  called  the  time  grid.  Both  lower  and  upper  bounds  on  the  lengths  of  the 
periods  are  assumed.  We  note  that  the  existence  of  a  finite  time  grid  is  usually  assumed  in 
discrete-time  planning  and  control.  For  event-based  flows  included  in  the  model,  the  set  A 
indicates  the  possible  times  of  events. 

It  will  be  convenient  when  representing  inventory  balance  equations  to  work  with  cumu¬ 
lative  flows.  For  a  flow  x.  the  corresponding  cumulative  flow  is  denoted  by  the  upper  case 

t 

letter  V.  and  is  defined  by  X(t )  =  ^x(r)dr.  t  >0.  Conversely,  the  unintegrated  flow 

corresponding  to  the  cumulative  flow  X  is  denoted  by  the  lower  case  letter  x. 

Unless  otherwise  stated,  if  a  component  of  a  vector  of  flows  is  not  explicitly  defined, 
then  it  is  assumed  to  be  the  zero  function. 

Intermediate  Product  Transfers 

We  index  the  different  inputs  and  outputs  of  the  system  from  1  to  K  with  the  first  \1 
indices  denoting  products  or  materials  and  the  last  K  -  M  indices  denoting  non-storable  ser¬ 
vices  such  as  labor  and  machine  time.  Let  vt)  =  (v,j.vj.  ...  vff)  denote  the  vector  of  flows 

representing  transfers  of  products  from  activity  i,  /  =0,1 N,  sent  to  activity  j, 

j  =  1.2 . V  +  l.  These  transfers  may  not  be  immediately  received  at  activity  j.  To  distin¬ 

guish  transfers  sent  from  those  received,  the  framework  associates  with  vtJ  the  vector 

v,j  =  ( vj . v")  of  transfers  from  activity  /  received  at  activity  j .  For  example,  if  there  is 

a  constant  time  lag  ltJ  for  shipments  from  i  to  j .  then 

v,,U  -  /<,)  if  t  >  l,j 
prespecifted  otherwise. 

We  remark  that  the  transformation  between  transfers  sent  and  transfers  received  is  not 
modeled  as  a  production  activity  of  the  network.  If  such  a  transformation  requires  scarce 
resources  such  as  material  handling  equipment,  then  the  activity  network  should  be  re-defined 


to  represent  the  material  handling  as  an  activity. 

Allocations,  Transfers  and  Applications 

The  framework  further  distinguishes  transfers  received  at  an  activity  from  inputs  actu¬ 
ally  applied  in  production.  Let  u,  =  .w*1  *2 . u*)  denote  the  allocation  flows  of  sys¬ 
tem  exogenous  services  to  activity  /  and  let  y,  =  (y/.y,’ . y,K)  denote  the  flows  of  inputs 

actually  applied  in  production  at  activity  i. 

In  each  model,  a  joint  domain  D  for  the  allocations,  applications  and  transfers  ( ( w, } . 
(y, )  and  {v,y })  reflects  the  domain  of  technologically  feasible  flows.  The  domain  may  involve 
complicated  constraints  linking  applications  at  different  activities.  For  example,  suppose 
activities  .-I  and  B  represent  the  production  of  products  .4  and  B,  respectively.  Suppose  each 
activity  uses  the  service  of  the  same  machine  which  can  only  process  one  type  of  product  at  a 
time  Then  the  domain  for  applications  must  be  defined  to  ensure  that  applications  of  the 
machine  service  at  activities  .4  and  B  are  not  simultaneous. 

Dynamic  Production  Functions 

For  each  activity  /,  a  dynamic  production  function  /,  maps  input  flows 
.V,  =  (y,1.  .  .>,*)  applied  at  activity  i  into  realized  output  flows 

/,(.»•,)  =  (//(}',) . /,w(y,  ))•  The  domain  of  the  production  function  f,  (as  derived  from  the 

joint  domain  D)  is  denoted  by  D,.  Typically,  one  or  more  components  of  y,  will  be  the  zero 
function.  Note  that  we  depart  from  the  traditional  economic  definition  of  a  production  func¬ 
tion  in  two  respects.  First,  our  production  functions  map  applications  into  realized  outputs 
rather  than  allocations  into  possible  outputs  Second,  our  production  functions  do  not  relate 
steady-state  rates  of  inputs  and  outputs,  but  rather  they  map  functions  of  time  representing 
technologically  feasible  input  applications  into  functions  of  time  representing  realized  outputs. 

(  he  Model  of  Production 

Given  an  identification  of  the  model  elements  and  their  domains,  the  framework 


specifies  a  deterministic  model  ot  production  in  terms  ot  llow  ^insi-ruimn  .tiM'j  ■  i 
lows.  Let  //"(/)  denote  the  inventors  at  time  t  at  aciivitv  <  ot  produ^i  i  ".j!-  <  j  ■> 

1</m<A/.  At  each  time  r  >  0.  the  cumulative  suppls  of  product  or  material  <>•  jij,  •. 

including  production,  received  transfers  and  initial  stock  -  is 

V 

O0)  -  /•lm( >, ►</ )  -  2,*7,n-  ■-  ■  - 

;  "0 

and  the  cumulative  amount  of  product  or  material  m  transferred  from  and  used  bv  actmtv  . 
is 

,V.| 

2  *  >/"</).  <2  21 
j-t 

For  each  product  or  material  m.  inventors  balance  constraints  require  that  for  each  activity  ; 
and  all  time  t  >  0. 

.v  v  ♦  i 

i?(n  =  /,m<0)  «•  /••,'"(>•,)(/)  +  -  2  >,7(o  -  Y,mtn  *  o.  <2. 3) 

)' 0  1  * 1 

We  refer  to  (2.3)  as  the  fundamental  inventory  balance  equation  The  framework  models 
inventories  of  products  at  each  activity  allowing  for  completed  products  awaiting  transfer  and 
transferred  products  awaiting  application.  Hence,  a  particular  product  may  have  multiple 
inventory  locations.  It  is  important  to  note  that  material  in  process,  either  in  transit  between 
activities  or  inside  an  activity  production  process,  is  not  considered  to  be  inventory. 

On  the  other  hand,  for  each  non-storable  input  m.  \l  +  \  <m  <K ,  we  require  for  each 
activity  t  and  all  time  t  >0  that 

yru)<»rv).  c.4) 

Equation  (2.4)  expresses  the  conservation  of  service  m  by  activitv  i. 

Let  u  =  Ut'.u1.  .  uu)  denote  a  vector  of  final  output  flows.  For  each  product  m  and 

all  time  t  >0.  inventory  balance  constraints  at  the  sink  require  that 

,v 

/.v . i (0)  ^  21*«"!v.i<o  -  L'mv)  -  °- 


(2.5) 
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Finally,  the  framework  models  conservation  of  exogenous  inputs  as  follows.  Let  cm(t) 
denote  the  system  supply  or  capacity  of  exogenous  input  m  at  time  t .  If  input  m  is  storable, 
then  the  cumulative  allocation  may  not  exceed  the  cumulative  supply,  i.e.,  for  all  time  t  >0, 

,v 

2l'S|,(0<Cm(/).  (2.6) 

i  =  i 

If  input  m  is  non-storable,  then  for  all  time  t  >0. 

2  <"(/)<cm(/).  (2.7) 

i  =  i 

In  the  case  of  dynamic  processes,  activity  output  flows  and  transfers-received  flows  from 
time  0  until  some  time  after  time  0  could  be  functions  of  applications  and  transfers-sent 
which  occur  at  or  before  time  0.  We  assume  that  all  flows  which  are  consequences  of  applica¬ 
tions  or  transfers-sent  at  or  before  time  0  are  fixed  and  prespecified  in  planning  calculations. 

From  the  perspective  of  the  general  framework,  a  deterministic  model  of  a  specific  pro¬ 
duction  process  is  categorized  and  explicitly  described  by  identifying  and  defining  the  abstract 
model  elements:  the  activity  network,  the  exogenous  inputs,  the  intermediate  and  final  pro¬ 
ducts,  the  activity  dynamic  production  functions,  and  the  domains  for  the  allocations,  appli¬ 
cations  and  transfers.  While  the  activity  network,  inputs  and  products  are  well-defined  in 
most  presentations  of  the  models  we  consider  in  follow-on  sections,  the  production  functions 
and  the  domains  are  typically  only  implicitly  treated.  We  will  thus  focus  our  attention  on 
these  model  elements. 

Rather  than  repeatedly  delineate  assumptions  common  to  various  models,  it  is  con¬ 
venient  to  define  the  following  model  categories.  In  an  acyclic  production  model,  the  activity 
network  is  acyclic.  In  a  product-generated  production  model,  each  activity  produces  a  single 
product  not  produced  by  any  other  activity.  In  a  discrete-time  production  model,  each  rate- 
based  flow  is  a  step-function  constant  in  each  time  period.  In  addition,  a  finite-horizon  is 
assumed.  In  a  normal  production  model,  (l)  activities  do  not  receive  products  they  can  not 
apply,  and  (2)  activities  do  not  transfer  products  they  can  not  produce. 
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3.  LINEAR  PROGRAMMING  MODELS 

We  now  illustrate  applications  of  the  framework  for  developing  more  accurate  models  of 
production.  Starting  with  explicit  sets  of  assumptions  about  the  physical  process,  we  use  the 
framework  to  derive  linear  programming  models  of  production.  We  show  that  familiar  formu¬ 
lations  incorporating  time  lags  do  not  correctly  model  the  physical  phenomena  they  purport 
to  represent.  We  provide  a  correct  formulation,  and  we  also  generalize  the  formulation  to 
accommodate  non-integer  lead  times  and  unequal  length  time  periods,  cases  which  have  not 
been  treated  before. 

3.1.  Leontief  Production  Functions 

Most  familiar  linear  programming  models  of  production  employ  a  very  simple  class  of 
activity  production  functions.  In  particular,  if  the  applications  of  inputs  are  positive,  then 
they  must  be  proportional  and  therefore  may  be  indexed  in  terms  of  one  profile.  The  outputs 
produced  are  also  assumed  proportional  and  indexed  by  the  same  profile  which  indexes  the 
inputs.  We  refer  to  the  profile  as  the  intensity  of  the  activ  ity. 

Casting  this  in  the  language  of  the  framework,  we  say  that  the  domain  D,  of  a  produc¬ 
tion  function  /,  is  Li  ontief  if  there  are  non-negative  constants  a,m.  m  -  1.2 . K  such  that 

D,  c  { v,  |  y,mU )  =  a,m:,(t),  m  =  1.2 . A.'  for  some  intensity  curve  } 

As  each  curve  y,m  in  a  Leontief  domain  D,  may  be  expressed  in  terms  of  its  defining  intensity 
curve  we  write  instead  of  /,(>*,).  We  say  that  a  production  function  /,  with  Leontief 

domain  is  itself  Leontief  if  there  are  non-negative  constants  c,m,  m  =  1.2, . \t  such  that 

/,"■(;,)  =  c,mr, ,  m  -  1,2 . \f .  (The  production  functions  are  so  named  due  to  their  simi¬ 

larity  to  the  steady-state  functions  used  in  Leontief  [1951].)  The  coefficient  u,"1  represents  the 
rate  input  m  is  applied  to  A,  per  unit  intensity,  and  the  coefficient  e”  represents  the  rate  pro¬ 
duct  m  is  produced  by  I,  per  unit  intensity. 

If  an  activity’s  production  process  is  modeled  by  a  Leontief  production  function,  then 
the  model  assumes  production  is  instantaneous  and  time-invariant.  That  is,  an  output  rate  of 
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an  activ  ity  at  a  point  in  time  is  solely  a  function  of  input  application  rates  at  the  same  point 
in  time:  moreover,  this  input/output  transformation  is  the  same  at  all  points  in  time.  If  out¬ 
put  at  time  t  is  actually  a  function  of  input  applications  over  some  period  of  time,  then  the 
Leontief  function  is  inappropriate. 

3.2.  Basic  Linear  Programming  Models  Without  Time  Lags 

The  following  assumptions  about  the  physical  system  lead  to  a  linear  programming 
model  of  production: 

( 1)  The  system  is  represented  by  a  product-generated  (normal)  activity  network. 

(2)  Each  activity  dynamic  production  function  is  Leontief  and  the  only  constraints  on  its 

domain  D,  are  those  implied  by  the  Leontief  assumption.  (This  characterizes  the 
Dynamic  Linear  Activity  Analysis  Model  (DLL4.U)  presented  in  Shephard  et  al.  [1977].) 
As  a  result  of  assumption  (1),  activity  intensity  will  be  measured  in  units  of  output. 
Thus.  /,'(--)  =  -.  1-1.2 . M. 

(3)  There  are  no  domain  constraints  on  allocations  or  transfers  (other  than  non-negativity). 

(4)  A  product  is  held  in  inventory  only  at  the  activity  which  produces  it. 

(5)  Intermediate  product  transfers  are  "instantaneous’,  i.e..  transfers  sent  equal  transfers 
received.  Combining  this  assumption  with  assumptions  (3)  and  (4),  as  intermediate  pro¬ 
ducts  are  transferred,  they  are  applied.  In  the  language  of  the  framework, 

v1  ssv*  =  i* *  *  d 1  z 

-  vij  -  aj^j- 

<6)  Transfers  of  final  product  are  instantaneous,  fn  view  of  assumptions  (3)  and  (4), 
transfers  of  final  product  are  equated  to  final  output  requirements,  i.e..  vj Ni.t  =  u‘ . 

(7)  The  only  exogenous  inputs  are  the  non-storable  services,  i.e.,  vo!=0,  m  =  \ . W. 

/  =  1 . .V. 

(8)  Since  domain  constraints  on  allocations  are  assumed  to  be  trivial,  allocations  are 

equated  to  application,  i.e..  w*  =  >*.  k  =  At  + 1 . K. 
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From  assumptions  ( 1  )-(8).  we  obtain  the  following  model  of  production:  For  each  intermedi¬ 
ate  product  i,  inventory  balance  equation  (2.3)  reduces  to 

St 

/,'(/)  =  //( 0)  +  Z,U)  -  Zoj'ZjU)  -  L'U)  >  o.  all  t  >0.  (3.1) 

j- i 

t 

For  each  non-storable  service  k,  the  allocation  constraint  (2.7)  becomes 

s 

^  ck{t).  all/>0.  (3.2) 

I  »  1 

To  develop  a  linear  programming  model,  we  invoice  the  additional  assumptions  of  a 
discrete-time  production  model  with  unit  length  time  periods  (0,1],  (1,2],  .  (As  is  cus¬ 

tomary,  we  refer  to  the  interval  (/-!,/]  as  "period  t~.)  Note  that  on  each  unit  interval,  all 
flows  in  and  out  of  inventory  of  product  /  are  at  constant  rates;  hence  the  rate  of  change  of 
the  inventory  level  during  each  interval  is  constant.  If  we  enforce  nonnegativity  of  the  inven¬ 
tory  level  of  product  i  at  the  end  points  of  an  interval,  it  follows  that  the  inventory  must  be 
nonnegative  during  the  entire  interval.  Hence,  inventory  balance  can  be  guaranteed  for  all 
t  >0  simply  by  enforcing  (3.1)  at  the  time  grid  points. 

Equations  (3.1)  and  (3.2)  now  may  be  rewritten  as.  respectively. 

I  i  st 

l‘(t)  =  //( 0)  +  2  "«»  -  S  (  +  ^  o,  1-1.2 . W,  f  =  1.2,  ,  (3.3) 

r  -  I  t  -  I  J  - 1 

St 

2a,*-„  <  cfa.  k  =  Xf  + 1 . K,  t  =  1,2,  •  •  ,  (3.4) 
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amount  of  service  k  required  per  unit  of  product  i,  and  aj  denotes  the  amount  of  product  / 


input  per  unit  output  of  product  j.  We  have  just  shown  that  physical  assumptions  ( 1 )-( 8 ) 
lead  to  the  constraints  (3.5)-(3.7)  in  the  traditional  linear  programming  multi-period  multi¬ 
stage  aggregate  production  planning  model  (see,  for  example,  Hax  and  Candea  [1985]  or 
Johnson  and  Montgomery  [1974]). 

The  assumptions  under  which  (3.5K3.7)  comprise  a  valid  model  of  production  are  now 
clear.  If  the  transformations  of  activity  input  to  output  are  not  proportional,  not  time- 
invariant,  or  not  instantaneous,  or  if  transfers  become  significant  either  because  there  are 
shipment  lags  or  because  transfers  are  not  applied  immediately  (e.g.,  inspections),  then  (3.5)- 
(3.7)  is  no  longer  a  valid  model  of  production.  (We  remark  that  assumptions  (2),  (3),  (5)-(8) 
are  necessarily  assumed  when  one  posits  the  model  (3.5)-(3.7);  however,  alternatives  to  (1) 
and  (4)  exist  which  also  lead  to  (3.5M3.7).) 

3.3.  Linear  Programming  Models  With  Time  Lags 

One  particularly  unrealistic  aspect  of  the  assumptions  leading  to  the  previous  linear  pro¬ 
gramming  model  is  the  requirement  that  processes  occur  instantaneously.  For  example, 
assumption  (2)  insists  that  output  from  an  activity  depends  only  on  input  to  that  activity  at 
the  same  instant  in  time.  This  assumption  precludes  operations  with  significant  post¬ 
processing  lags  (e.g.,  time  for  steel  to  cool),  and  also  operations  with  significant  processing 
times  (e.g.,  metal-cutting).  Likewise,  assumption  (5)  precludes  situations  with  significant 
transfer  times. 

Considerable  research  has  been  devoted  to  developing  practical  linear  programming 
models  without  these  instantaneousness  assumptions.  In  all  of  this  work,  all  significant  time 
lags  are  lumped  under  the  single  heading  'lead-time'.  For  example,  Billington  el  al.  [1986] 
describe  lead  time  as  a  'nonproduction  lag,  such  as  the  time  for  paint  to  dry.  hot  metal  to 
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I  Ml  I 

/'(/)  =  /'( 0)  +  2-.(r -L.)  -  2  2 ajzj,  -  2«r  ^  0,  /  =  1,2,  .  (3.8) 

r» l  j - 1 r» I  r» 1 

Here.  L,,  the  "unavoidable  lead  time"  for  product  /.  is  defined  so  that  output  of  product  i 
started  in  period  t  may  first  become  available  in  period  t  +  L,.  Inherent  in  the  definition  of 
L,  is  the  assumption  that  each  unavoidable  lead  time  is  independent  of  the  production  rate 
and  can  be  expressed  as  an  integral  number  of  time  periods. 

We  now  demonstrate  that  these  parameters  can  not  correctly  model  the  stated 
phenomena.  To  do  so.  we  define  three  different  types  of  lags  contributing  to  lead  time  and 
formulate  a  correct  linear  programming  model  assuming  these  time  lags  are  integral. 


Output  Lag 

In  some  production  processes,  a  product  can  not  be  transferred  or  released  to  inventory 
for  a  period  of  time  after  it  is  "produced".  For  example  after  painting,  a  fixed  amount  of  time 
may  be  required  for  the  paint  to  dry;  or  after  lumber  is  sawed,  a  fixed  amount  of  time  may  be 
required  to  inspect  and  grade  the  output.  For  such  processes,  assumption  (2)  must  be 
modified  as  follows;  Input  applications  at  -1,  are  still  proportional,  but  the  output  curve  is 
shifted  by  a  constant  amount  L4 ,  >0.  (The  symbol  ‘LA’  denotes  "lag  after"  application).  In 
the  language  of  the  framework. 

:,(t  -  L- 1, )  if  r  >  L.4, 

prespecified  for  0 <t<LA, 

Note  that  applications  of  inputs  are  still  instantaneous,  but  consequent  output  emerges  after  a 
time  lag  LA,.  Note  also  that,  since  applications  are  prespecified  before  time  0.  output  of  pro¬ 
duct  i  must  be  prespecified  for  re(0, LA,].  Substituting 

i 

F',(Z,)U)  =  |/,'(-Hr)</r  =  |  Z,(r)dr. 


//(-,  Hr)  = 
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Transfer  Lag 

A  period  of  time  may  be  required  to  transfer  parts  between  activities.  Suppose  a 
transfer  from  activity  /  to  activity  j  requires  exactly  LTtJ  time  units.  (The  symbol  ‘LT’ 
denotes  "lag-transfer  .)  Retaining  our  other  assumptions,  assumption  (5)  must  be  modified  as 
follows: 

v!j(t  ~  LTtJ)  if  t>LT,j 

\j(t)  -  dj-j(r)  -  i,j(r)  -  peespecified  for  0<T^£.ry. 

Note  that  z,(t)  must  be  prespecified  for  re(0.L7"y]  because  it  corresponds  to  transfers  sent 
from  A,  at  or  before  time  0.  Retaining  the  output  lags  and  noting  that 


—  -  i/ 

J  *'(, 


(2.3)  now  becomes 


»  <+lt„ 


/,'</)  =  //( 0)  +  /  :,(T)dr  -  2  /  -  i‘U)  >  0,  all  t  >0.  (3. 

-Li,  i- 1  tr„ 

The  intensity  of  each  .1*  must  now  be  prespecified  for  -Hk  <t  <max{Z.7",*  }. 


Input  Lag 

Another  type  of  lag  arises  when  there  is  a  delay  between  the  time  material  is  received  at 
an  activity  and  the  time  it  is  applied.  For  example,  there  may  be  an  inspection  of  inputs 
required  which  does  not  utilize  scarce  resources.  Alternatively,  such  a  lag  may  be  used  to 
approximate  queue  time  required  because  of  cyclic  batch  production  on  a  single  machine,  or 
other  complex  domain  constraints.  We  further  modify  assumption  (5)  so  that  there  is  a  lag  of 
exactly  LB,,  >0  time  units  between  transfer  receipt  at  A,  of  intermediate  product  /  and  its 
application  at  Ar  (The  svmbol  ‘LB'  denotes  "lag  before"  application.)  In  the  language  of  the 
framework  the  modified  assumption  (5)  is 

[vy(r  -  LB,j)  =  vy(r  -  LT,,  -  LB:I )  if  r>  LTtJ  +  LBI) 


\j{r)  -  (jj.jir)  | prespecified 


for  0  <  t  <  LT,j  ♦  LB,j 


Retaining  our  other  modifications  and  observing  that 


t+LT,+LB,, 


17,(0  =  /  ‘V,(r)r/r  =  j  aj:j(T)dT, 

LT,  LT„*LB„ 


(2.3)  now  becomes 

t  -  LA  t  +LT  *  LB  , 

/,'(/)  =  //«))  +  /  r,(r)Jr  -  2  {  aj:j(T)dT  -  L'U)  >  0,  all  t  >0.  (3.11) 

With  the  additional  lag,  the  intensity  of  each  .4*  must  now  be  prespecified  for 
-L4*  <  r  <  max  { LT,/,  +  LB,i<  )  . 

I 

It  is  important  to  note  that  the  fixed  lag  LB,,  between  transfer  receipt  and  application 
predetermines  a  quantity  a):,(T)  of  inventory  of  product  /  residing  at  A,  during  the  interval 
(t-LB,j,t\.  Hence.  //(/)  no  longer  represents  the  total  system  inventory  of  product  /  at  time 
t,  but  rather  the  portion  of  total  inventory  which  is  i.ot  stock  committed  by  the  lag  assump¬ 
tions.  (We  do  not  bother  to  write  the  inventory  balance  constraint  for  Ij(t)  since  it  is 
automatically  satisfied  by  our  assumptions.) 

To  develop  a  linear  programming  model,  we  once  again  invoke  the  assumptions  of  a 
discrete-time  production  model  with  unit-length  time  periods.  If  Li,  and  (LTtJ  +  LB,,)  are 
integers,  then  all  flows  in  and  out  of  inventory  of  product  i  are  constant  during  each  interval. 
Recalling  our  earlier  discussion,  inventory  balance  will  be  guaranteed  for  all  />0  if  it  is 
enforced  at  the  time  grid  points.  Retaining  the  notation  :,T  for  the  constant  intensity  of  A, 
and  u'j  for  the  constant  rate  of  demand  for  product  i  during  (r-  l,r],  the  discrete-time  version 
of  (3.1 1)  is 

I  Ml  1 

l‘,(l)  =  /,'( 0)  +  -  S  2) a):i{r*LT ,*LB., )  -  -  0-  1  =  1.2.  •  (3.12) 

r ■ l  / -  1  r»l  r -  1 


Inspecting  the  form  of  (3.11)  and  (3.12).  it  appears  that  separate  identification  of 
transfer  LTt)  and  input  LB,,  lags  is  unnecessary.  This  is  true  under  the  special  assumptions 
we  have  made  in  this  model:  (1)  the  only  domain  constraints  on  transfers  are  nonnegativity. 


and  (2)  uncommitted  inventory  of  a  product  is  held  only  where  produced.  Under  alternatives 
to  (I).  such  as  discrete  domains  for  transfers  or  upper  bounds  on  transfers,  it  is  meaningful  to 
relax  (2).  In  such  a  case,  uncommitted  inventory  of  an  intermediate  product  would  be  held 
both  where  produced  and  where  used,  necessitating  separate  transfer  and  input  tags  as  well  as 
multiple  inventory  balance  equations  to  ensure  conservation  of  the  product. 

Comparison  with  Familiar  Formulations 

We  contrast  (3.12)  with  familiar  formulations  incorporating  time  lags.  The  inventory 
balance  constraints  in  the  mathematical  programming  model  of  Billington  et  al.  [1983,  1986] 
are  given  in  (3.8).  Hax  and  Candea  [1985]  also  use  (3.8)  without  providing  any  physical 
definition  of  "production  lead  time".  Billington  et  al.  [1983,  1986]  claim  their  lead  time  L, 
models  such  phenomena  as  time  to  transfer  parts  or  time  for  paint  to  dry.'  However,  any 
time  to  transfer  parts  from  intermediate  product  inventory  to  an  activity  which  will  apply  the 
product  as  input  is  a  lag  of  the  form  LTtJ.  It  is  not  correct  to  include  allowances  for  lags  of 

the  form  LT,t  or  LBt)  in  the  authors’  parameter  L,  in  (3.8).  a  point  we  now  elaborate. 

Consider  the  simple  example  in  Figure  1.  Application  of  a  quantity  of  product  i  at 
activity  j  corresponds  to  resource  application  at  activity  z  to  produce  that  quantity  at  least  6 
time  periods  earlier.  Similarly,  application  of  product  i  at  activity  k  corresponds  to  resource 
application  at  activity  z  at  least  10  time  periods  earlier.  If  T,  =  6,  then  there  are  choices  for 
the  :„'s  and  r*/s  feasible  in  (3.8)  but  physically  impossible.  If  L,  =  10,  then  there  are  choices 
for  the  :,fs  and  z/s  physically  possible  but  infeasible  in  (3.8).  If  6  <  L,  <  10,  then  both 
problems  arise.  On  the  other  hand,  there  can  be  no  doubt  if  one  uses  (3.12)  with  LA,  =  3, 

LT,j  +  LB,j  =  3,  LT +  LB,k  =  7,  since  (3.12)  properly  defines  the  physically  possible 

choices  for  the  :„'s.  z/s.  and  zk,'s . 

Even  if  it  were  the  case  that  LT,t  +  LBtJ  =  LTlk  +  LBlk.  (3.8)  would  still  be  invalid. 
The  issue  in  this  case  is  precisely  what  the  variables  l'(t)  and  the  parameter  /'( 0)  in  (3.8) 
represent.  Although  not  explicitly  defined  in  the  papers  cited  above,  a  common  interpretation 
of  /'(/)  in  (3.8)  is  that  it  represents  the  total  system  quantity"  of  product  t  which  is  "in 
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Figure  1 

Example  Network  with  Lead  Times. 
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advance  of  use,"  i.e.,  in  advance  of  application  at  another  activity  or  withdrawal  from  the  sys¬ 
tem  by  external  demands.  From  our  definition  of  the  elements,  this  interpretation  of  /'(/)  is 
actually  given  by 

l  Ml  t 

ru)  =  /'< o)  +  2-*  <r  -LAt)-  2  2 a‘j=AT)  -  2w^- 

r*  1  >  =  I  r*  I  t*  1 

which  corresponds  to  (3.8)  if  one  takes  L,  =  Li,.  However,  if  one  includes  a  positive 
allowance  for  LT,j  +  LB,t  in  L ,,  then  l‘(t )  in  (3.8)  no  longer  represents  system  inventory  of 
product  z  at  time  /,  nor  for  that  matter  any  other  measure  of  inventory  of  product  i  at  time 
t .  It  would  be  impossible  to  properly  represent  initial  inventory  conditions  in  (3.8). 

In  general,  any  model  with  transfer  or  input  lags  which  uses  an  initial  "total  system 
inventory"  parameter  /'  (0)  is  flawed:  initial  conditions  are  not  correctly  expressed.  That  is, 
/'( 0)  will  in  general  include  amounts  of  product  z  to  be  immediately  applied  at  a  follow-on 
activity  Ar  amounts  still  being  inspected,  amounts  in  transit,  and  amounts  not  yet  transferred 
from  .-I,.  When  there  are  shipment  or  inspection  lags,  the  total  quantity  /'( 0)  should  not  be 
simultaneously  made  available  to  follow-on  activities.  On  the  other  hand,  in  (3.12)  the 
parameter  //( 0)  properly  represents  only  inventory  of  product  z  not  yet  committed  to  follow- 
on  activities,  while  the  parameters  of  the  form  cyl,  :j2.  .  .  .  .  serve  to  specify  the 

status  of  quantities  of  product  z  already  sent  to  activity  .4,  which  are  in-transit,  being 
inspected,  and  to  be  immediately  applied.  1 

1  The  correct  expression  for  the  total  system  quantity  /'(f)  of  product  i  at  time  z  is  derived  as  follows. 

/'(z)  is  defined  by 

U 

/'(z)  -  //</ )  +  2i/;(z)*/^(z)i, 

i  • i 

where  ITij(t)  denotes  the  quantity  of  product  z  in-transit  between  A,  and  A,  at  time  z.  In  the  language  of 
the  framework. 


j 

i 


i;u)  -  /;< o)  +  f>)  -  y;u).  and  /77,(z)  *  rr„( o>  ♦  »v,(z>  -  v;,u) 

Our  lag  assumptions  imply  that 

/;( 0)  -  |  a;z;(r)dr  and  IT]{ 0)  -  |  ajz^T  . 

Upon  substitution  and  simplification  /'(z)  is  expressed  as 

/'(/)-//( />♦  2i  2  .,1.  z-i.2. 
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Krajewski  and  Ritzman  [1977]  propose  a  "general  model"  which  utilizes  balance  equa¬ 
tions  of  the  form 


w 

/'(/)  -  l'(t  -  1)  =  -  “l  (3.13) 

j- 1 

which  is  a  different  special  case  of  (3.11).  In  this  case.  L,  cannot  reasonably  include  an 
allowance  for  phenomena  such  as  time  for  paint  to  dry  on  product  i.  Any  allowance  for 
(Li,  +  LT,j  +  LB,j)  included  in  L,  must  be  independent  of  i.  If  in  fact  (LA,  +  LT,,  +  LB,,) 
varies  with  i.  then  their  model  incorrectly  restricts  the  intensity  of  .A,  relative  to  the  intensi¬ 
ties  of  its  immediate  predecessor  activities.  Unlike  the  previous  special  case,  this  error  in  the 
restriction  of  intensities  can  be  overcome  by  replacing  the  authors’  parameter  L,  in  (3.13) 
with  parameters  LtJ  =  (LA,  +  LT,,  +  LB,j).  2 

However,  once  again  initial  conditions  can  not  be  expressed  correctly.  In  the  model 
defined  by  (3.13),  product  /  enters  inventory  immediately  following  resource  application  but 
before  paint  is  dry.  The  single  parameter  /' (0)  in  (3.13)  can  not  account  for  the  age  distribu¬ 
tion  of  inventory  (e  g.,  dry  parts  versus  those  still  drying).  On  the  other  hand,  in  (3.12)  the 
parameter  /,'( 0)  properly  includes  only  dry  parts  ready  for  transfer  to  follow-on  uses,  while  the 
parameters  -h2-la  )-  •  -koi  specify  the  status  of  parts  which  are  still  drying. 

Our  final  comment  of  the  comparison  with  familiar  formulations  concerns  the  notion  of 
a  "frozen  planning  horizon"  (Hax  and  Candea  [  1 985]).  This  notion  recognizes  that  produc¬ 
tion  rates  over  a  certain  initial  period  of  time  must  be  prespecified  in  planning  calculations 
Relative  to  the  measurement  of  initial  inventory  at  time  0,  Billington  cl  al.  [1983.1986]  and 
Hax  and  Candea  [1985]  prespecify  production  intensities  of  .-1,  during  the  interval  (-f.,.0). 
while  Krajewski  and  Ritzman  [1977]  prespecify  production  intensities  of  A,  during  the  inter¬ 
val  (0 ,Lj],  The  former  approach  is  not  correct  unless  L,  =  L.A,  and  LT,,  +  LBtJ  =  0  for  all  j\ 
the  latter  approach  is  not  correct  unless  L,  =  LT,,  *  LB,,  for  alt  i  and  LA,  =  0.  As  we  have 

seen,  in  general,  one  must  prespecify  the  intensity  of -1,  during  the  interval 
:  A  similar  improvement  may  be  made  to  MRP  calculations.  See  Section  4. 
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(-LA,,  ma \{LTJt  +  LB Jt }], 

as  well  as  measure  /,'( 0),  in  order  to  commence  planning  calculations  at  time  0. 

A  Linear  Programming  Model  With  Non-Integer  Lead  Times 

The  assumption  that  lead  times  are  integer  is  quite  restrictive,  although  it  is  assumed  by 
all  the  familiar  formulations.  We  now  derive  a  discrete-time  model  when  lead  times  are  not 
necessarily  integral.  We  retain  the  assumption  that  activity  intensities,  service  capacities,  and 
final  demands  are  constant  rates  during  unit  length  time  periods. 

To  evaluate  inventory  balance  of  product  /  at  time  t .  one  must  express  the  integral  func¬ 
tions  in  (3.1  1)  in  terms  of  the  ,'s  and  :/r's  in  a  manner  which  accounts  for  the  non-integer 
limits  of  integration.  For  a  real  number  v,  we  let  v ’  denote  the  smallest  integer  greater  than 
or  equal  to  v.  and  let  \  ~  denote  the  largest  integer  less  than  or  equal  to  x.  The  first  integral  in 
(31  |  >  is  expressed,  for  all  /  >  0.  as 


r 


U):‘<-ur  if t -LA,  <(-L4,y 

[( -LA,  r  +  H,  ]rl(  ,L<  +  2  -<r  (3.14) 

i  -n  r  <rsi(  -Li  i  otherwise. 

-  [<(  -LA,)-U  -LA,V\ L^- 


In  (3.14),  the  coefficients  of  the  first  and  last  terms  of  the  lower  expression  express  the  frac¬ 
tions  of  intensity  in  the  first  and  last  time  periods  which  are  included  within  the  limits  of 
integration;  the  middle  term  simply  sums  up  intensities  of  all  time  periods  (if  any)  in 
between.  The  upper  expression  accounts  for  the  degenerate  case  when  both  limits  of  integra¬ 
tion  lie  in  the  same  time  period.  Next,  let  LX,t  =  LTtj  +■  LB,r  The  second  integral  in  (3.1 1)  is 
analogously  expressed,  for  all  /  >0.  as 


■IX 


(rk/r 


r  if  t+LX.j  <(LX„r 

[LX,;  -LX„\z,lLX *  2  -V  (3.15) 

u.x  r  <fsn -lx.  >  otherwise. 

*■  [((  +LX,j)-U  +  LX,j)~]:j(i  t  LX  y 


Finally.  C  '<()  in  (3.1  1 )  is  expressed,  for  all  t  >  0.  as 
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«d 


m 


i  '(/)  =  +  (r  -  /-)«/-. 


(3.16) 


Substituting  the  identities  (3.  l4)-(3. 16)  into  (3.11)  provides  the  desired  expression  of  //(f)  in 


terms  of  discrete-time  variables  for  any  /  SO.  Identifying  the  Hows  which  are  functions  of 


applications  or  transfers-sent  occurring  before  time  0.  we  note  that  for  each  .1*  and  for 


(-/..(*)*.( -/_!*)*  «■  I.  .0.1. 


max  ( LA'*  }. 


the  intensities  (rtr)  must  be  prespecified. 


If  the  time  lags  are  not  integer,  ensuring  inventory  balance  (3.11)  only  at  the  integer 


time  grid  points  is  not  enough  to  ensure  feasibility.  Consider  the  example  shown  in  Figure  2. 


Here.  L  ‘  =0,  LA,  =  1.7.  /.  V,(  =  1.7.  and  </,'  =  !.  A  simple  check  shows  that  (3.1 1)  is  satisfied 


for  t  =  1.2.3.  However,  at  t  =  1.7  cumulative  output  of  /  is  zero,  yet  two  units  are  required  by 


j  at  this  time--  a  clear  infeasibility.  This  phenomenon  can  happen  since  the  flows  //(-, ),  v,.. 


and  vtJ  are  no  longer  constant  on  the  given  time  intervals.  These  flows  are  still  step  functions. 


but  the  time  points  where  rates  of  flow  may  change  fall  in  between  the  integer  time  points. 


Ma.iy  investigators  overcome  this  difficulty  by  either  overlooking  it  and  rounding  lead 


times  to  integer  amounts  or  by  subdividing  the  natural  lime  period  so  that  lead  times  become 


integral.  The  first  approach  either  under-estimates  lead  time  -  potentially  leading  to  infeasi¬ 


bilities  -  or  over-estimates  lead  times  -  potentially  leading  to  excessive  work-in-process.  The 


second  approach  greatly  increases  problem  size. 


To  ensure  inventory  balance  (3.1 1)  at  all  time,  we  claim  it  is  necessary  and  sufficient  to 


check  (3.1  I)  at  all  time  points  when  the  rates  of  final  demand  flows,  activity  output  flows  or 


intermediate  product  transfer  flows  may  change  For  product  i  let  T,  denote  this  collection  of 


time  points.  A  point  in  time  t  >  0  is  an  element  of  f,  if.  and  only  if.  t  -  LA,  or  t  *■  LX,,  or  r 


itself  is  an  integer.  In  the  example. 


T,  =  (0.  3.  7.  10.  13.  1.7.  2  0.  2.3,  2.7.  3.0. 


For  proof  of  the  claim,  note  that  on  the  intervals  defined  by  successive  points  of  T,.  all  flows 


'VV  .  •/? 


.>'■  >  ■/*  o'.V, 
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in  and  out  of  inventory  of  product  i  are  at  constant  rates;  hence,  following  our  earlier  discus¬ 
sion.  the  inventory  level  at  any  time  point  within  such  an  interval  will  be  nonnegative  if  the 
inventory  level  is  enforced  to  be  nonnegative  at  the  end  points. 

This  formulation  is  not  only  accurate  but  it  requires  significantly  fewer  constraints  than 
would  be  included  in  a  model  employing  a  time  grid  fine  enough  to  make  all  lead  times 
integer.  In  the  example,  our  formulation  requires  3  balance  equations  per  period  for  product 
i.  To  realize  integer  lead  times  a  grid  with  periods  of  length  0.10  is  required  for  product  /  , 
making  for  7  extra  balance  equations  within  each  period.  Worse,  if  one  time  grid  is  used  to 
enforce  inventory  balance  for  all  products,  an  even  finer  grid  might  be  required.  The  width  of 

SI 

the  required  grid  is  precisely  the  least  common  divisor  of  the  points  in  UTj.  On  the  other 

i-i 

hand,  the  precision  of  the  lead  times  does  not  affect  the  number  of  balance  equations 
required  in  our  proposed  approach. 

A  Linear  Programming  Model  With  Unequal  Length  Time  Periods 

Another  restrictive  assumption  of  familiar  formulations  is  that  all  time  periods  have 
equal  length.  However,  unequal  length  periods  are  desirable,  for  example,  when  natural  time 
periods  such  as  weeks  include  varying  numbers  of  working  days  because  of  vacations  or  holi¬ 
days.  We  now  modify  our  previous  formulation  to  handle  the  case  when  the  lengths  of  the 
periods  with  constant  activity  intensities  are  not  necessarily  equal.  Let  S  denote  the  set  of 
epochs  marking  the  end  points  of  these  periods.  To  ensure  feasibility,  it  is  necessary  and 
sufficient  to  evaluate  (3.1 1)  at  those  time  points  t  such  that  t  -  L4,  or  t  +  L\,t  or  /  is  an  ele¬ 
ment  of  5  (thus  modifying  the  definition  of  T,).  The  expressions  of  the  integrals  in  (3.14), 
(3.15)  and  (3.16)  under  these  general  conditions  must  be  modified  since  the  coefficients  of 
intensity  variables  and  demand  rates  were  derived  assuming  a  period  length  of  one.  To 
obtain  the  correct  discrete-time  version  of  (3.11),  we  modify  (3.14),  (3.15)  and  (3.16)  as  fol¬ 
lows.  First,  for  a  real  number  x,  re-define  x*  (resp.  x~)  to  mean  the  smallest  (resp.  largest) 
epoch  tk  eS  not  less  than  (resp.  not  greater  than)  x.  Next,  substitute  the  ratios 
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t  l-L4,)*  +  LA.  (t-LAg)-  (t -L4j)~ 

- — — — — - - — ,  -  and  -  (3.17) 

1,)-  (~L4,r  ~(-L4,y  (t-L4,)* -U-LAiY 

for  the  coefficients  (/).  [(-£.4,)*  +  LA,]  and  [(/ -L4,)  -  (t -L.4,)']  in  (3.14),  respectively.  (We 

define  -jj  =0.)  The  denominators  of  the  ratios  in  (3.17)  express  the  lengths  of  the  time  periods 

in  which  the  limits  of  integration  lie;  the  ratios  thus  express  the  fraction  of  intensity  in  those 
periods  which  lies  within  the  limits  of  integration.  Likewise,  we  substitute  the  ratios 


t  LX*  -  LXtJ  and  U  +  LX^-V+LXjj)- 

LX.J-LX,-'  LX* -LX-  (t+LXtJ)*  -{t  +  LX.j)- 


(3.18) 


for  the  coefficients  (r).  [LX,J  -  LX,j]  and  [(t  +LXt))  -  (t  +LX,j)~]  in  (3.15),  respectively. 
Finally,  we  modify  (3.16)  to  become,  for  all  /  >0. 


I'M/) 


+ 


(/-/-) 
(/+  -  /■) 


(3.19) 


Incorporating  the  modifications  (3.17),  (3.18)  and  (3.19),  inventory  balance  constraints  (3.11) 

enforced  at  the  set  of  epochs  T,  for  product  i,  /  =  1 . V/,  and  capacity  constraints  (3.4) 

enforced  for  each  epoch  in  5  form  a  more  general  linear  programming  production  model.  The 

M 

oserall  time  grid  for  this  discrete-time  model  is  A  =  U  T,. 

i-i 

In  sum,  we  have  developed  an  accurate  yet  parsimonious  linear  programming  model  for 
product-generated  networks  with  three  different  types  of  lead  times  for  each  product.  A  time 
grid  (with  possibly  unequal-length  periods)  for  resource  allocation  is  assumed  to  be  given. 
Activity  input  applications  are  constant  on  these  intervals  according  to  Leontief  domains. 
Based  on  the  lead  times,  distinct  sets  of  epochs  for  each  product  are  established  for  enforcing 
inventory  balance. 


Our  more  general  formulation  is  by  no  means  the  most  general  linear  programming  pro¬ 
duction  model.  The  lead  times  in  the  foregoing  model  serve  to  relax  the  assumption  of  instan¬ 
taneousness  of  intermediate  product  applications,  outputs,  and  transfers  which  is  assumed  by 
linear  programming  models  without  time  lags.  The  assumption  of  instantaneous  applications 


of  serv  ice  resources  can  also  be  relaxed,  i.e.,  the  assumption  of  Leontief  domains  is  not  neces¬ 
sary  for  obtaining  a  linear  programming  model.  For  example,  Krajewski  and  Ritzman  [1977] 
propose  a  formulation  with  resource  input  coefficients  which  distribute  resource  application 
over  an  integral  number  of  time  periods  between  intermediate  product  application  and  reali¬ 
zation  of  output.  Using  the  framework,  we  could  extend  their  model  to  admit  noninteger 
durations  for  resource  applications. 

As  a  final  remark,  if  there  are  significant  processing  times  (i.e.,  noninstantaneous  service 
applications),  these  times  also  must  be  accounted  for  in  the  lags  separating  intermediate  pro¬ 
duct  input  applications  and  consequent  output  of  each  activity.  It  is  not  correct  to  delete  pro¬ 
cessing  times  from  the  lags  appearing  in  inventory  balance  equations,  as  has  been  proposed  by 
Billington  et  al.  [1986].  Moreover,  setup  times  introduce  more  breakpoints  in  inventory  input 
and  output  flows;  if  setup  times  are  fractional,  enforcing  inventory  balance  only  at  integer 
times  does  not  guarantee  inventory  balance  throughout  continuous  time.  A  full  development 
of  linear  programming  models  for  the  case  of  activity  production  functions  with  noninstan¬ 
taneous  service  applications  is  beyond  the  scope  of  this  paper;  however,  linear  programming 
models  for  particular  noninstantaneous  cases  of  interest  are  developed  in  Hackman  and 
Leachman  [1985a,  1985b,  1986]  and  Leachman  [  1986], 


4.  MANUFACTURING  RESOURCES  PLANNING 


Manufacturing  Resources  Planning  (MRP)  systems  convert  time-phased  final  output 
requirements  into  requirements  for  intermediate  products  and  raw  materials,  and  into 
requirements  for  service  resources.  The  logic  of  these  conversions  is  described  as  a  "time- 
phased  explosion"  through  a  "bill  of  materials'  structure  (Orlicky  [1975]).  In  this  section,  we 
redefine  MRP  in  terms  of  the  framework.  MRP  is  seen  as  a  schedule  of  intermediate  product 
transfers  rather  than  a  schedule  of  production.  We  provide  a  simple  analysis  of  why  the  use 
of  MRP  may  result  in  excessive  work-in-process  inventories.  Alternatives  and  modifications 
to  standard  MRP  logic  are  discussed. 

The  MRP  Model 

Consider  an  acyclic,  product-generated,  normal,  discrete-time  production  model  with 
unit-length  time  periods.  There  are  no  exogenous  materials  in  the  model.  (All  raw  materials 
of  interest  are  included  as  products  in  the  network.)  The  activities  (products)  are  ordered 
such  that  if  .4,  supplies  intermediate  product  to  ,4y,  then  i  <j.  Final  transfer  requirements 
l',„v+|.  (  =  1.2 . iV.  are  specified  for  the  system.  (Since  the  model  is  product-generated,  we 

N*l 

suppress  the  superscripts  on  the  flows.)  Let  V,(t)=  2  KjU)  denote  the  cumulative  total 

y-.  +  l 

transfer  requirement  of  .1,  at  time  t.  Let  R,(t)  =  max  (0,F,(/)  -  /'(0)}.  In  MRP  parlance, 
v,(f)  is  the  "gross  requirement"  for  product  i  at  time  l,  while  r,{t)  is  the  "net  requirement"  for 
product  i  at  time  /.  Final  transfer  requirements  ( vl  V4.,(f))  are  termed  the  "master  production 
schedule".  In  MRP  systems,  gross  and  net  requirements,  and  the  master  production  schedule 
are  expressed  as  event-based  flows  on  an  equal-length  discrete  time  grid  such  as  weeks. 

MRP  systems  include  a  material  requirements  planning  (mrp)  module  to  determine 
requirements  for  products  and  materials.  Coefficients  [a])  are  prespecified,  indicating  the 
amount  of  product  i  required  as  input  to  produce  one  unit  of  product  j.  So-called  "lead 
times'  (L, }  are  also  prespecified.  Each  L,  is  an  integer  which  defines  a  time  offset  between 
the  event  when  a  transfer  requirement  for  product  i  is  due  and  the  event  when  intermediate 
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product  inputs  used  to  produce  the  requirement  are  withdrawn  from  inventory  by  A,.  For 

t  =  V..V-  I . t  and  for  all  /.  the  following  computations  are  performed  recursively: 

.v,i 

K(l)  =  2  KlU)  (4.1) 

/ »l  ♦  I 

R,(i)  =  max  {0.r<(/)-/,(0)}  (4.2) 

=  a/R,(t  +  L, ).  all  j  <i.  (4.3) 

An  additional  step  is  sometimes  included  between  (4.2)  and  (4.3)  to  modify  the  flow  R,(t)  to 
reflect  desirable  lot  sizes.  The  (modified)  flows  ( r,  )  are  the  result  of  mrp  calculations;  each 
flow  r ,  defines  order  quantities  of  product  i  due  at  the  time  grid  points. 

MRP  systems  also  include  a  capacity  requirements  planning  (crp)  module  to  determine 
"requirements"  for  service  resources.  Coefficients  { zz*>  are  prespecified,  indicating  the  amount 

of  service  k  required  as  input  to  produce  one  unit  of  product  i,  i  *  1,2 . N.  Following 

execution  of  the  mrp  module,  the  crp  module  estimates  the  system  requirement  for  service  k 
in  period  t  as 

,v 

(4.4) 

i-i 

If  (4.4)  exceeds  capacity  of  service  k  in  period  f .  then  one  is  obliged  to  reduce  in  some  way 
the  final  transfer  requirements  v,  v „ , ,  /  =  1,2 . .V. 

Analysis  of  the  mrp  Module 

MRP  systems  incorporate  an  incomplete  model  of  production.  Activity  application 
flows,  activity  production  functions  and  domain  constraints  are  not  defined.  In  general,  the 
requirements  curve  R ,  for  activity  t(  is  run  the  cumulative  output  curve  F,(Y,)  of  .4,  simi¬ 
larly,  the  transfer  curve  V)t  is  not  the  cumulative  curve  Y{  of  application  of  product  /  at  .4,. 
In  order  for  MRP  to  be  valid,  fundamental  inventory  balance  at  each  activity  .4,  must  be 
ensured,  i.e., 

F,{Y,)>R,  (4.5) 

and 


Instead  of  specifying  V,  and  F,  ()',).  MRP  systems  replace  (4.5)  and  (4.6)  with  the  relationship 
(4.3).  Attempting  to  ensure  feasibility  solely  using  (4.3)  can  lead  to  excessive  work-in-process 
inventory,  as  we  now  discuss. 

For  the  purposes  of  illustration,  we  consider  application  of  MRP  to  a  production  system 
satisfying  the  assumptions  of  Section  3  2.  We  further  assume  that  the  rate  of  production  (per 
period)  by  activity  .3,  is  bounded  by  For  example.  5,  may  reflect  an  allocation  of  service 
resources  to  .4,.  Plotted  in  Figure  3  is  a  particular  net  requirements  curve  R,.  Also  plotted  is 
the  latest"  output  curve  Z,L  satisfying  the  requirements  R,.  Z,L  is  generated  by  assuming  it  is 
feasible  to  produce  up  to  5,  in  each  period.  We  plot  the  case  in  which  r,U)  exceeds  5,  in  vari¬ 
ous  periods,  as  is  frequently  the  case  in  actual  practice. 

Suppose  one  were  to  produce  item  i  according  to  the  "late"  schedule  dictated  by  the  out¬ 
put  curve  Z,L  shown  in  Figure  3.  Because  the  intermediate  products  necessary  for  production 
must  be  available  prior  to  production,  if  item  /  is  required  in  production  of  item  i.  then 
I),  >  Y/  =  a/  Z,L.  To  minimize  the  inventory  of  item  j  at  4,  one  would  take  I  ;|  =  a{ Z,L. 
Note  that  one  must  calculate  Zf*  from  R,  in  order  to  obtain  the  l';1  which  minimizes  inven¬ 
tory.  (A  simple  algorithm  for  calculating  Z,L  from  R,  in  a  discrete-time  model  is  given  in 
Leachman  (1979].) 

However,  in  an  MRP  system.  V,,  is  calculated  using  (4.3).  Ideally,  the  time  shift  used 
should  be  the  smallest  such  shift  guaianteeing  that 

i)1(/)><i//1tin.  (4  7) 

In  this  example.  L,  =  4  is  the  minimal  such  time  shift.  (See  Figure  3.)  In  general,  the  smallest 
L,  is  dependent  on  the  particular  net  transfer  requirements  curve  R,  ,  and  the  bound  on  the 
rate  of  production  f, .  As  discussed  before,  [.,  is  exogenously  prespecified,  independent  of  R , 
and  5,  A  large  L,  must  be  chosen  to  ensure  feasibility.  Even  if  were  the  mimmum  shift 
satisfying  (4  7).  the  insistence  on  a  simple  time  shift  model  as  expressed  in  (4  3)  in  general 


Time 


3 
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results  in  excessive  intermediate  product  inventories,  as  depicted  in  Figure  4. 


In  Figure  4.  a{  =  I.  and  the  cumulative  transfer  of  product  j  to  A,  is  taken  to  be 
R,(t  +  4)  (i.e.,  a  time  shift  of  required  output).  Suppose  production  of  item  i  proceeds  at  the 
maximum  rate  allowable  by  the  bound  on  the  production  rate  and  the  rate  of  intermediate 
product  /  transferred  to  it.  The  resulting  cumulative  output  curve  is  labeled  Z,£  (the  "early 
curve  ).  By  following  Z,E  instead  of  Z,L.  the  inventory  of  intermediate  product  j  is 
significantly  reduced.  However,  the  inventory  of  item  /  awaiting  transfer  becomes  significant. 
Actual  production  must  lie  somewhere  in  between  the  ZE  and  ZE  curves.  Hence,  when  a  sim¬ 
ple  time  shift  model  is  used  to  relate  transfers  in  lieu  of  production  functions,  inventory  of 
product  awaiting  production  or  inventory  of  product  awaiting  transfer  will  in  general  be  large. 

To  reduce  excess  inventory,  one  must  calculate  output  curves  from  requirements  curves, 
and  then  calculate  transfers  from  output  curves.  If  service  resources  are  preallocated  to 
activities  (as  in  the  example),  these  calculations  are  easy  to  incorporate  into  the  mrp  calcula¬ 
tion  (Leachman  [1979]).  Alternatively,  one  can  try  to  solve  the  mathematical  programming 
problem  representing  simultaneous  resource  allocation  and  scheduling  (Billington  et  al 
[1983.  1986]). 

As  discussed  in  Billington  el  al.  (1983,  1986],  the  programming  formulation  is  computa¬ 
tionally  intractable  for  realistic  problems.  The  mrp  calculation  remains  a  popular  method  for 
large-scale  scheduling;  moreover,  inflation  of  the  lead  time  parameters  is  a  convenient  means 
of  addressing  stochastic  aspects  of  production  which  are  difficult  to  model.  Nonetheless,  a 
simple  improvement  to  the  mrp  calculation  can  be  made  which  is  analogous  to  the  improve¬ 
ment  we  suggest  for  l.p.  formulations  incorporating  lags.  Consider  application  of  MRP  to  a 
production  system  satisfying  the  assumptions  of  Section  3.3.  The  mrp  parameter  L,  must 
account  for  time  to  transfer  inputs  from  predecessors  ( LTn ),  time  to  inspect  inputs  ( LB ,,). 
time  to  produce  (LA,),  as  well  as  time  spent  in  inventory  /,'  (i.e..  time  resulting  from 
differences  between  output  and  requirements  curves,  as  discussed  above).  The  first  two  fixed 
lags  depend  on  the  source  of  inputs;  hence  the  mrp  parameter  L,  in  (4.3)  must  include  an 


allowance  for  ma \(LT n  +  LB„)  in  addition  to  other  factors.  In  the  case  of  multiple  activity 

inputs  with  different  lags,  some  reduction  in  lead  time  (and  excess  inventory)  is  obtained  by 
substituting  for  L,  in  (4.3)  a  parameter  LJt  which  is  calculated  the  same  as  L,  except  that  the 
allowance  above  is  replaced  by  (LT^  +  LBJt ). 

Analysis  of  the  crp  Module 

We  now  turn  to  the  representation  of  service  resources  in  MRP.  From  the  point  of  view 
of  the  framework  it  is  not  valid  to  estimate  service  resource  loads  in  terms  of  transfer  require¬ 
ments,  as  is  done  by  the  crp  module.  As  we  have  seen,  even  when  activity  production  func¬ 
tions  are  Leontief  the  requirements  curve  for  an  activity  A,  is  not  its  intensity  curve. 

An  infeasibility  suggested  by  the  crp  module  does  not  necessarily  imply  that  the  final 
transfer  requirements  are  infeasible,  nor  even  that  the  MRP-derived  transfer  schedules  are 
infeasible.  Consider  once  again  our  example  in  Figure  4.  The  dashed  curve  is  the  continuous 
flow  corresponding  to  the  event-based  R,(t)  curve.  The  crp  module  estimates  resource  con¬ 
sumption  by  .4,  in  terms  of  the  dashed  curve.  Since  its  slope  exceeds  z,  (the  slope  of  Zf-(t))  in 
several  periods,  the  crp  module  would  indicate  that  R,(t)  is  infeasible;  however.  Figure  4 
clearly  shows  that  R,(t)  is  feasible  with  respect  to  5,. 

We  have  seen  that  MRP  systems  incorporate  a  model  of  production  in  which  all  flows 
are  expressed  in  terms  of  intermediate  product  transfers.  On  the  other  hand,  basic  linear  pro¬ 
gramming  models  incorporate  a  model  of  production  in  which  all  flows  are  expressed  in  terms 
of  intensities  of  activity  input  application.  In  a  more  accurate  model  of  multi-stage  produc¬ 
tion,  it  is  likely  that  both  application  and  transfer  flows  must  be  explicit  in  the  model  in  order 
to  represent  the  various  domain  constraints  peculiar  to  each  type  of  flow. 
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5.  CRITICAL  PATH  MODELS 

In  familiar  critical  path  models  (CPM).  an  acyclic  network  of  activities  .4  { . 4  v  is 

given.  Each  .-I,  is  assumed  to  operate  in  some  un-interrupted  interval  of  time  with  integer 
length  d,  >0.  In  the  activity-on-node  format,  an  arc  from  A,  to  A,  means  there  is  a  strict  pre¬ 
cedence  relationship  between  .1,  and  Aj.  (,4;  can  not  be  started  until  after  A,  has  finished.)  In 
the  standard  application  of  resource-constrained  CPM,  resource  use  by  an  activity  is  assumed 
constant  during  its  duration  (Moder  et  al.  [1985]). 

Resource-constrained  CPM  embeds  a  model  of  production  as  follows.  The  activity  net¬ 
work  is  the  CPM  activity-on-node  network,  exclusive  of  the  source  for  exogenous  inputs  and 
the  sink  for  final  outputs.  The  underly  ing  production  system  has  the  following  characteristics. 

(1)  The  production  model  is  acyclic  and  normal  with  a  finite  horizon. 

(2)  Each  activity  produces  a  vector  of  outputs,  whereby  a  distinct  intermediate  product  is 
supplied  to  each  follower.  Activities  with  no  followers  produce  a  final  product 
transferred  to  a  sink  node.  Since  we  may  identify  uniquely  the  products  of  the  system 
with  arcs  in  the  network,  we  replace  the  notation  m  for  product  with  the  notation  (/,;). 

(3)  Exactly  one  unit  of  each  product  is  produced. 

(4)  Transfers  are  event-based  flows.  Quantity  transferred  equals  0  or  1.  We  take  the 
time  grid  .V  to  be  the  set  of  nonnegative  integers. 

(5)  Activity  production  functions  have  a  restricted  Leontief  domain.  Input  applications  to 
.4,  are  indexed  by  a  "box'  curve  of  the  form 

d,  if  t6(5,,S(  +  d,\ 

0  otherwise 

where  S,  e  A  corresponds  to  the  start  time  of  A,.  Domain  constraints  on  allocations  are 
assumed  to  be  trivial. 

(6)  Cumulative  output  of  an  activity  is  measured  in  terms  of  the  fraction  of  required 


resources  which  have  been  applied.  We  define  =  for  each  j  such  that  (i,j) 

corresponds  to  an  arc  in  the  network.  Note  that  an  activity  simultaneously  produces 
one  product  for  each  follow-on  activity. 


It  is  a  simple  matter  to  verify  that  strict  precedence  is  implied  by  the  fundamental 
inventory  balance  equation  and  the  assumptions  (l)-(6).  Hence  strict  precedence  between 
activities  is  simply  a  particular  form  of  inventory  balance  which  arises  when  there  are  event- 
based  transfers. 

The  representation  (l)-(6)  of  CPM  in  terms  of  the  framework  changes  according  to  one’s 
conventions  for  describing  the  physical  phenomena.  Assumption  (6)  is  based  on  a  particular 
convention  for  output  measurement,  e  g.,  when  50%  of  the  required  effort  has  been  made,  an 
activity  is  considered  50%  done.  A  second  representation  arises  if  activity  output  is  viewed  as 
a  discrete  lump"  emerging  exactly  when  resource  applications  are  complete.  In  this  represen¬ 
tation.  we  relax  the  integer  constraint  on  vftJ)  in  (4)  and  replace  (6)  by 


f,u  j\ )(t) 


I  if  r  =  S,  +  d,  and  (/.))  is  an  arc 
0  otherwise. 


(5.2) 


Now.  each  production  function  has  Leontief  domain  but  is  no  longer  Leontief. 

A  third  representation  arises  if  instead  of  the  output  or  transfers  being  discrete,  the 
follow-on  applications  of  intermediate  product  are  viewed  as  discrete.  In  this  case,  we  main¬ 
tain  the  original  assumption  (6),  relax  the  domain  for  transfers  in  (4)  to  simple  nonnegativity, 
and  alter  (5)  such  that 


y},J)(r ) 


1  if  r  =  Sj  and  (i.j )  is  an  arc 
0  otherwise. 


(5.3) 


Now,  the  domains  for  the  production  functions  are  no  longer  Leontief:  only  service  resource 
applications  (>•*)  are  indexed  by  an  intensity  function  (5.1). 

In  fact,  a  strict  precedence  model  may  be  formulated  without  event-based  flows.  In  this 
fourth  model,  there  is  no  inventory  of  product  U.j)  at  .4,,  i.e.,  (4)  becomes  v,<jj)  =  f,{,J\y,). 
Each  production  function  has  restricted  Leontief  domain  exactly  as  in  (5).  With  these 


assumptions,  we  can  maintain  strict  precedence  if  we  replace  (6)  by 


Hr) 


dt  if  r  e  <.S,  +•  d,  .S,  +  d ,  +  dj]  and  ( t  .j )  is  an  arc 
0  otherwise. 


(5.4) 


The  production  function  defined  in  (5.4)  maps  the  intensity  curve  for  .4,  into  the  earliest 
intensity  curve  for  consistent  with  the  assumption  of  strict  precedence.  That  is.  the  produc¬ 
tion  function  determines  the  appropriate  bound  on  the  choices  for  intensity  of  Aj.  This  last 
representation  is  the  only  one  of  the  four  which  is  preserved  under  aggregation  of  resource- 
constrained  CPM  networks;  using  this  model,  a  computationally  tractable  approach  to  multi- 
project  aggregate  planning  has  been  developed  (Hackman  and  Leachman  [1985a,  1985b], 
Leachman  and  Boysen  [1985]).  See  Figure  5  for  a  pictorial  representation  of  the  four  alterna¬ 
tive  models. 

Comparing  CPM  to  the  basic  l.p.  model  in  Section  3.2.  note  that  both  models  utilize  the 
same  production  function,  except  that  application  flows  in  CPM  have  the  severely  restricted 
domain  defined  by  (5.1).  It  is  natural  to  consider  relaxation  of  these  domain  constraints  to 
allow  flexibility  in  activity  operation  akin  to  that  allowed  by  the  activity  analysis  model  in 
linear  programs.  In  many  project-oriented  physical  systems,  the  application  of  service 
resources  need  not  be  at  one  fixed  rate  from  start  to  finish  of  project  activ  ities.  (See,  for  exam¬ 
ple,  the  discussion  of  the  application  of  trade  labor  services  to  ship  overhaul  activities  in 
Leachman  [1979]  and  Leachman  [1983].) 

A  number  of  authors  have  developed  extensions  in  this  regard,  which  we  now  categorize 
in  terms  of  the  allowed  domains  for  application  flows.  Wiest  [1967]  and  Talbot  [1982]  have 
considered  resource-constrained  scheduling  when  there  are  discrete  alternatives  for  the  dura¬ 
tion  (and  associated  resource  requirements)  of  each  activity.  However,  from  the  point  of  view 
of  the  framework,  each  alternative  still  corresponds  to  a  'box’  intensity  curve  for  each 
activity,  i.e..  a  constant  rate  of  resource  application  from  start  to  finish  of  an  activity.  Leach¬ 
man  [1983]  develops  a  technique  for  resource  leveling  when  there  is  a  continuous  range  of 


(1)  Leontief  production  function,  event-based  transfer 
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possible  intensity  levels  for  each  activity,  i.e.,  (5.1)  is  replaced  by 

-4-  if  re  (5,.  5,  +  d,\  for  some  S,  e  R  »  and  -j-  e  [r. .5, ) 

U|  dt 

"  0  otherwise. 

Weglarz  [1981]  and  Dincerler  [1984]  relax  this  domain  even  further  in  considering  resource- 
constrained  scheduling.  They  develop  scheduling  algorithms  which  exploit  the  assumption 
that  intensity  may  vary  during  the  duration  of  an  activity,  i.e..  (5.1)  is  replaced  by 

:  for  some  5,.  F,  e  R+,  z,  is  nonzero  only  on  (5.6) 

the  interval  [S, .F, ],  and  for  re [5, ,F,  ],  z,  <  z,{t)  <  z,)  . 

In  summary,  we  have  cast  CPM  as  a  model  of  production,  recognizing  that  each  arc  in 
the  CPM  network  represents  the  transfer  of  a  distinct  intermediate  product,  and  that  a  strict 
precedence  relationship  is  an  expression  of  inventory  balance  of  such  a  product.  Our  simple 
categorization  of  extensions  to  CPM  in  terms  of  relaxations  of  the  domain  for  applications 
illustrates  the  power  of  the  framework  to  compare  models  easily  yet  rigorously. 


6.  A  STRUCTURED  APPROACH  TO  MODELING  PRODUCTION 

We  have  introduced  a  general  framework  for  deterministic  models  of  production  and 
illustrated  its  use  through  the  evaluation  and  improvement  of  familiar  models.  However,  the 
chief  intent  of  our  framework  is  to  guide  the  formulation  of  specific  models  of  production. 
By  “specific'  we  mean  that  a  particular  physical  production  process  is  considered,  the 
planning/scheduling  decision  levels  are  specified,  and  one  seeks  to  construct  a  model  of  pro¬ 
duction  for  each  level. 

The  process  of  model-building  using  the  framework  proceeds  in  an  orderly,  structured 
fashion  as  follows: 

Step  1 :  Develop  a  detailed  model  of  production. 

(a)  Elucidate  assumptions  about  the  physical  system. 

(b)  Identify  the  model  elements  (production  functions,  domain  constraints,  etc.). 

(c)  Express  constraints  which  ensure  conservation  through  continuous  time  of  all  services, 
materials,  and  products  of  interest. 

(d)  Approximate  as  necessary. 

Step  2:  Develop  an  aggregate  model  of  production  for  the  decision  level  under  consideration 
(if  required). 

Repeat  steps  (a)-(d)  above  for  the  aggregate  level  of  detail. 

Step  3:  Test  and  validate  the  model. 

Note  that  we  develop  an  aggregate  model  appropriate  for  the  decision  requirements  of 
the  planning  problem  having  first  developed  an  accurate  model  of  the  production  system.  We 
view  aggregation  as  a  meta-model  operation  rather  than  a  model-specific  operation.  In  our 
framework,  aggregation  requires  extensive  modeling  effort.  It  leads  from  one  (detailed)  model 
consistent  with  the  general  framework  to  another  (aggregate)  model  consistent  with  the  frame¬ 
work.  In  our  experience,  these  two  models  may  -  and  often  do  -  have  dramatically  different 
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structure.  See.  for  example,  Leachman  and  Boysen  [1985]  and  Hackman  and  Leachman 
[1985a.  1985b].  On  the  other  hand,  previous  treatments  of  aggregation  consider  only  model- 
specific  aggregation  which  is  a  largely  mathematical  operation  for  reducing  the  dimension  of  a 
model  while  preserving  its  structure.  It  leads  from  one  instance  of  a  model  (the  detailed- 
level)  to  another  instance  (the  aggregate-level)  having  the  same  structure.  For  example,  row 
and  column  aggregation  begins  with  a  linear  program  whose  primitive  elements  (.4,  b,  c.  and 
v)  are  already  specified.  It  leads  -  without  additional  modeling  -  to  another,  albeit  smaller, 
linear  program.  (We  do  not  mean  to  say  that  model-specific  aggregation  techniques  are  not 
valuable.  In  fact,  they  may  be  useful  components  of  solution  techniques  for  aggregate  models 
developed  via  our  structured  approach.) 

In  this  paper  we  have  illustrated  only  Step  1  of  the  structured  approach,  by  way  of 
analysis  of  the  familiar  models.  We  derived  linear  programming,  MRP  and  CPM  models 
using  the  framework  from  explicit  assumptions  about  the  physical  system.  Examples  of  the 
development  of  novel  aggregate  production  models  using  the  structured  approach  and  the 
framework  may  be  found  in  Hackman  and  Leachman  [1985a,  1985b.  1986]. 
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